Comparison of multigrid and incomplete LU shifted-Laplace preconditioners for the inhomogeneous Helmholtz equation∗
نویسندگان
چکیده
Within the framework of shifted-Laplace preconditioners [Erlangga, Vuik, Oosterlee, Appl. Numer. Math., 50(2004), pp.409–425] for the Helmholtz equation, different methods for the approximation of the inverse of a complex-valued Helmholtz operator are discussed. The performance of the preconditioner for Helmholtz problems at high wavenumbers in heterogeneous media is evaluated. Comparison with other preconditioners from the literature is also presented.
منابع مشابه
Domain decomposition preconditioning for high-frequency Helmholtz problems with absorption
In this paper we give new results on domain decomposition preconditioners for GMRES when computing piecewise-linear finite-element approximations of the Helmholtz equation −∆u − (k2 + iε)u = f , with absorption parameter ε ∈ R. Multigrid approximations of this equation with ε 6= 0 are commonly used as preconditioners for the pure Helmholtz case (ε = 0). However a rigorous theory for such (so-ca...
متن کاملShifted-Laplacian Preconditioners for Heterogeneous Helmholtz Problems
We present an iterative solution method for the discrete high wavenumber Helmholtz equation. The basic idea of the solution method, already presented in [18], is to develop a preconditioner which is based on a Helmholtz operator with a complex-valued shift, for a Krylov subspace iterative method. The preconditioner, which can be seen as a strongly damped wave equation in Fourier space, can be a...
متن کاملGMRES-based multigrid for the complex scaled preconditoner for the indefinite Helmholtz equation
Multigrid preconditioners and solvers for the indefinite Helmholtz equation suffer from non-stability of the stationary smoothers due to the indefinite spectrum of the operator. In this paper we explore GMRES as a replacement for the stationary smoothers of the standard multigrid method. This results in a robust and efficient solver for a complex shifted or stretched Helmholtz problem that can ...
متن کاملOn a class of preconditioners for solving the Helmholtz equation
In 1983, a preconditioner was proposed [J. Comput. Phys. 49 (1983) 443] based on the Laplace operator for solving the discrete Helmholtz equation efficiently with CGNR. The preconditioner is especially effective for low wavenumber cases where the linear system is slightly indefinite. Laird [Preconditioned iterative solution of the 2D Helmholtz equation, First Year’s Report, St. Hugh’s College, ...
متن کامل3D Helmholtz Krylov Solver Preconditioned by a Shifted Laplace Multigrid Method on Multi-GPUs
We are focusing on an iterative solver for the three-dimensional Helmholtz equation on multi-GPU using CUDA (Compute Unified Device Architecture). The Helmholtz equation discretized by a second order finite difference scheme is solved with Bi-CGSTAB preconditioned by a shifted Laplace multigrid method. Two multi-GPU approaches are considered: data parallelism and split of the algorithm. Their i...
متن کامل